Interaction between static holes in a quantum dimer 
model on the kagome lattice 

Gregoire Misguich, Didina Serban and Vincent Pasquier 

Service de Physique Thcorique, CEA Saclay, 91191 Gif-sur-Yvette Cedex, FRANCE. 
URA 2306 of CNRS 

Abstract. 

A quantum dimer model (QDM) on the kagome lattice with an extensive ground- 
state entropy was recently introduced [Phys. Rev. B 67, 214413 (2003)]. The ground- 
state energy of this QDM in presence of one and two static holes is investigated by 
means of exact diagonalizations on lattices containing up to 144 kagome sites. The 
interaction energy between the holes (at distances up to 7 lattice spacings) is evaluated 
and the results show no indication of confinement at large hole separations. 
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1. Introduction 

Quantum dimer models (QDM) can provide some effective descriptions of the low-energy 
singlet dynamics of frustrated quantum antiferromagnets T]. A basis of the Hilbert 
space of these models is made by all (nearest-neighbor) dimer coverings of the lattice 
and the Hamiltonian allows these dimer to move along local resonance loops. Two kinds 
of phases are well understood for these models in two dimensions: dimer crystals P and 
resonating valence-bond (RVB) dimer liquids [21 E] ■ Crystals are characterized by long- 
ranged dimer-dimer correlations and spontaneous lattice symmetry breaking and RVB 
liquids show no broken symmetry but topological order and Z2- vortex excitations |3]. 
Importantly also, these two phases can be distinguished by the behavior of holes (or 
spinous) when the system is "doped", that is we allow sites which are not occupied 
by any dimer. Holes experience a mutual interaction which grows linearly with their 
separation in a dimer crystal and this interaction confine them in pairs. On the other 
hand, they propagate independently in a RVB liquid background. A first and simple 
step is to consider QDM with static holes (non-magnetic impurities or spinous). In that 
case a relevant quantity is the ground-state energy as a function of the hole positions. 
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Table 1. The 8 different classes (up to rotations) of dimcrizations of an hexagon of 
the kagome lattice. 

This energy goes to a constant when two holes a separated far apart in a deconfined 
system whereas it grows hnearly in a confined system4 

We recently proposed a QDM on the kagome lattice which does not simply fall 
in any of these two categories 0. This model, hereafter called the /x-model, was 
introduced from the observation that the dimer kinetic energy terms arising from an 
overlap expansion of the spin-^ Heisenberg model PP generally have non trivial signs as 
soon as the competition of resonance loops with different lengths is considered jH] . The 
Hamiltonian is: 

n = -J2f^'^ (1) 

h 

where 

32 

IJ^h = Km +H.C (2) 

a=l 

I QDMs at Rokhsar-Kivelson points ^ are an exception: the ground-state energy remains exactly zero 
whatever the positions of the holes. 
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Figure 1. The hole forbids all resonance loops on hexagons 1 and 2 and suppresses 
some of the loops around 3 and 4. 



and h runs over the hexagons of the lattice and \da{h)) is one of the 32 possible 
dimerizations of h (table P). The sign counts the parity of the number of dimer 
involved§. It was realized that such signs can lead to a new state, different from dimer 
crystals or RVB liquids. In our previous study [S] the following results were obtained: i) 
The /i-model has an extensive ground-state entropy | log(2) per kagome site, that is 50% 
of the classical dimer entropy. This exponentially large degeneracy comes from a hidden, 
local, but non-abelian symmetry of the model, ii) It is possible to choose a basis of the 
ground-state manifold so that dimer-dimer correlations are short-ranged in each state. 
These ground-states are thus dimer liquids, iii) On the basis of exact diagonalizations 
we argued that, in addition to the ground-state degeneracy, the spectrum is likely to be 
gapless and that energy-energy correlations (as well as susceptibilities) are likely to be 
critical. 

In the present work we investigate numerically the effect of static holes in the /i- 
model. This issue is of particular importance as Dommange et al. [Hj pointed out in a 
recent work that static holes in the spin-i kagome antiferromagnet experience a short- 
distance repulsion and are probably deconfined at larger distances. Sindzingre et a/. [7] 
previously reached a similar conclusion about spinon deconfinement from an analysis of 
the value of the spin gap in a 24-site sample with two holes. We show in this paper that 
a somewhat similar behavior is observed in the /x-QDM. 

2. /i— model with holes 

As any QDM, the /i-model can be extended to include static holes. These holes can 
equally represent charge degrees of freedom or neutral spinous (unpaired spin in a dimer 
background). The new Hamiltonian Ti' contains all the kinetic terms of Ti except those 
where the resonance loop passes through a missing site. Consider a hole which belongs 
to two hexagons 1 and 2 (figure [l}. No loop of fii or /X2 survive in TC' because they would 
all pass through the missing site. As for hexagons 3 and 4, one half of their resonance 
loops pass through the hole and must be removed. In presence of a hole the operators 
/i3 and /i4 thus only contain 16 resonance loops (instead of 32). However these two 

§ In the absence of that sign the model reduces to that of reference [5] and can be solved exactly, it 
has a RVB liquid ground-state with topological order and gapped Z2-vortex excitations. 
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Figure 2. Energy cost A of a single hole in the fi QDM as a function of the number 
of hexagons Nh ■ 



modified operators satisfy tlie same algebraic relations IB/ as the hole-free 's. For any 
liexagon h 1,2 and for i = 3 or 4 we liave fif = 1 and: 

/Xj/i/i = fihfJ'i h not neiglibor of i (3) 

fJ'if^h = — f^hf^i h neighbor of i (4) 

These relations are easy to check with the help of the arrow representation of dimer 
coverings of the kagome lattice [01 El E|- It is also easy to check that the argument 
leading to an exponential degeneracy ~ 2^/® of the energy levels [5 holds even in the 
presence of these static holes. As a first result we thus find that the extensive ground- 
state entropy of the /i-model survives in the presence of holes. This also allows to 
use the reduced representation of the Hilbert space which was used in reference jlj to 
compute the spectrum in the absence of holes. The spectrum is non-degenerate in this 
representation, which has a dimension ~ 2^/^ (instead of 2^/^+^ for the dimer Hilbert 
space). The ground-state of systems up to 48 hexagons (144 kagome sites) can be 
obtained with a standard Lanczos algorithm. The result was checked (with and without 
holes) against direct calculations in the dimer Hilbert space for small systems (A^ < 48). 
We investigated samples with = 36, 48, 60, 72, 84, 108 and 144 kagome sites {Nh =12, 
16, 20, 24, 28, 36 and 48 hexagons). Periodic boundary conditions are used and the 
shapes of these clusters are the same as those of reference jH]. 

Interestingly this representation allows to compute the spectra of a even system 
pierced by a single hole. Strictly speaking the QDM is not defined on such an odd 
sample but the non-degenerate representation of the n algebra mentioned above can 
still be constructed. This trick is useful to estimate the energy cost A of a single hole 
in a given sample. 



3. Single hole ground-state energy 

The ground-state energy in the absence of holes is noted Eq, Ei is the energy with a single 
hole (two neighboring /i's removed). If {fJ^ifJ^j) correlations are neglected, the ground- 
state energy would increase by 2(/i) ~ 0.88 around each hole (the ground-state energy is 
estimated to be ~ —0.44 per hexagon in the thermodynamic limit In fact removing 
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Figure 3. Energy E2{d) of the ground-state of the /.t-QDM in presence of two holes 
at distance d. In the top panels this energy is compared to the energy Eq of the system 
without holes and in the lower panels E2{d) is compared to the hole-free energy E^y 
corrected by twice the energy cost A = £'i — £'o of a single hole. 



two neighboring jj, operators increases the ground-state energy by A = i^i — i?o ~ 0.6 
(figure |2|). It is easy to understand why the actual hole gap A is smaller than the 
naive estimate above. Because of the anti-commutation relations between nearby /i 
operators (equation 0} the system cannot simultaneously achieve a minimal energy [i.e 
/i = 1) on two neighboring hexagons. Removing some /i operators therefore decreases 
the frustration on their neighbors, which can acquire in turn a larger expectation value 
(lower energy). This larger "polarization" of the hexagons around the hole will enhance 
the frustration on their neighbors, and the corresponding /i will have to reduce (slightly) 
their expectation value compared with the bulk value. This mechanism produces spatial 
oscillations in (/i) (data not shown), oscillations which have the same wave- vector as 
the correlations which dominate in the bulk jS]. 

4. Two-holes ground-state energy 

The difference between the energy E2{d) with two holes at distance d and the energy 
Eq without holes is shown in figure El In the analysis of Ref. j3] it appeared that the 
/i-model has significant local (/ij/ij) correlations with a period of three hexagons. It 
is therefore convenient to plot separately the data for Nh not multiple of three and 
the others {Np = 12, 24, 36 and 48), which do not frustrate the local order. At short 
distance, when the two holes belong to a common hexagon, only 3 ^ operators are 
removed from 7i and the energy cost is roughly E2 — Eq ~ |A. This happens for c? = 1, 
-\/3 as well as for d = 2 when the two holes are on opposite sites of an hexagon. This 
is a short-distance effect because for d > 2 the number /i suppressed is always four. 
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In other words, when two holes sit on the same hexagon they minimize the number 
of loops which are "lost" for resonances ||. At intermediate distances (2 < c? < vT2) 
the energy decreases with distance in a regular way for all samples. In this range of 
d the behavior is thus reminiscent of the strong hole repulsion observed in the kagome 
Heisenberg model by Dommange et al. 0. 

For d> 4 the data suggest that the energy E2{d) goes to a constant. The values are 
indeed close to the energy Eo + 2A (see lower panels of figure Ej) which is expected if the 
dimer background was not mediating any interaction between the holes. We therefore 
argue that the /i model is not confining static holes. It is interesting to note that in 
the samples which do not frustrate the local order (right panels of figure E)) the hole- 
hole interaction seems to decay more slowly with distance than in the other samples. 
This rather slow decay is not incompatible with the interesting suggestion jHj of a l/d 
behavior. In addition, weak oscillations that can be observed. They may be related to 
the oscillations of (fi) mentioned in the previous section and to the (presumably) quasi 
long-ranged correlations in (fiifij)^ discussed in Ref. |S] 
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II This effect can also be found, to a smaller extent, in the two-hole energies of the spin-1/2 model |H]: 
Comparing the two ways two holes can be at distance d = 2, the energy is always lower when they 
belong to the same hexagon. There is, however, no strong reduction of E2 for d — 1 and \/3 in the spin 
model as we have in the /i-QDM. 



